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Abstract 



We classify Lie-Poisson brackets that are formed from Lie algebra extensions. The 
problem is relevant because many physical systems owe their Hamiltonian structure 
• to such brackets. A classification involves reducing all brackets to a set of normal 

forms, and is achieved partially through the use of Lie algebra cohomology. For 
extensions of order less than five, the number of normal forms is small and they 
' involve no free parameters. We derive a general method of finding Casimir invariants 

5^ ■ of Lie-Poisson bracket extensions. The Casimir invariants of all low-order brackets 

are explicitly computed. We treat in detail a four field model of compressible reduced 
magnetohydrodynamics. 
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1 Introduction 



This paper deals with the classification of Lie-Poisson brackets obtained from 
extensions of Lie algebras. A large class of finite- and infinite-dimensional dy- 
namical equations admit a Hamiltonian formulation using noncanonical brack- 
ets of the Lie-Poisson type. Finite-dimensional examples include the Euler 
equations for the rigid body [1] and the moment reduction of the Kida vor- 
tex [2] , while infinite-dimensional examples include the Vlasov equation [3,4] 
and the Euler equation for the ideal fluid [5-9]. Lie-Poisson brackets natu- 
rally define a Poisson structure (i.e., a symplectic structure) on the dual of 
a Lie algebra. For the rigid body, the Lie algebra is the one associated with 
the rotation group, 5*0(3), while for the Kida vortex moment reduction the 
underlying group is S0{2, 1). For the two-dimensional ideal fluid, the relevant 
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Lie algebra corresponds to the group of volume-preserving diffeomorphisms 
on the fluid domain. 

We will classify low-order bracket extensions and find their Casimir invariants. 
An extension is simply a new Lie bracket, derived from a base algebra (for 
example, S0{3)), and defined on n-tuples of that algebra. We are ruling out 
extensions where the brackets that appear are not of the same form as that of 
the base algebra. We are thus omitting some brackets [5,10], but the brackets 
we are considering are amenable to a general classification. 

The method of extension yields interesting and physically relevant algebras. 
Using this method we can describe finite-dimensional systems of several vari- 
ables and infinite-dimensional systems of several fields. For finite-dimensional 
systems an example is the two vector model of the heavy top [11]. For infinite- 
dimensional systems there are models with two [12-14], three [12,15,16], and 
four [17] fields. Knowing the bracket allows one to find the Casimir invari- 
ants of the system [18-20]. These are quantities which commute with every 
functional on the Poisson manifold, and thus are conserved by the dynam- 
ics for any Hamiltonian. They are useful for analyzing the constraints in the 
system [21] and for establishing stability criteria [22-26]. 

The outline of this paper is as follows. In Section 2, wc review the general 
theory behind Lic-Poisson brackets. We give examples of physical systems 
of Lie-Poisson type, both finite and infinite-dimensional. We introduce the 
concept of Lie algebra extensions and derive some of their basic properties. 
Section 3 is devoted to the more abstract treatment of extensions through 
the theory of Lie algebra cohomology [27-29]. Wc define some terminology 
and special extensions such as the semidirect sum and the Leibniz exten- 
sion. In Section 4, we use the cohomology techniques to treat the specific 
type of extension with which we are concerned, brackets over n-tuples. We 
give an explicit classification of low-order extensions. By classifying we mean 
reducing — through coordinate changes — all possible brackets to independent 
normal forms. We find that the normal forms are relatively few and involve 
no free parameters — at least for low-order extensions. In Section 5, we turn to 
the problem of finding the Casimir invariants of the brackets, those function- 
als that commute with every other functional in the algebra. We derive some 
general techniques for doing so that apply to extensions of any order. Some 
explicit examples are derived, including the Casimir invariants of a particular 
model of magnetohydrodynamics (MHD). These are also given a physical in- 
terpretation. A formula for the invariants of Leibniz extenions of any order is 
also derived. Then in Section 6 we use the classification of Section 4 to derive 
the Casimir invariants for low-order extensions. Finally in Section 7 we offer 
some concluding remarks and discuss future directions. 
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2 Lie— Poisson Brackets 



Lie-Poisson brackets define a natural Poisson structure on duals of Lie alge- 
bras. Physically, they often arise in the reduction of a system. For our pur- 
poses, a reduction is a mapping of the dynamical variables of a system to a 
smaller set of variables, such that the transformed Hamiltonian and bracket 
depend only on the smaller set of variables. (For a more detailed mathematical 
treatment, see for example [30-34].) The simplest example of a reduction is 
the case in which a cyclic variable is eliminated, but more generally a reduc- 
tion exists as a consequence of an underlying symmetry of the system. For 
instance, the Lie-Poisson bracket for the rigid body is obtained from a reduc- 
tion of the canonical Euler angle description using the rotational symmetry 
of the system [11]. The Euler equation for the two-dimensional ideal fluid is 
obtained from a reduction of the Lagrangian description of the fluid, which 
has a relabehng symmetry [26,35-37]. 

Here we shall take a more abstract viewpoint: we do not assume that the Lie- 
Poisson bracket is obtained from a reduction, though it is always possible to do 
so by the method of Clebsch variables [26] . Rather we proceed directly from a 
given Lie algebra to build a Lie-Poisson bracket. The choice of algebra can be 
guided by the symmetries of the system. After deriving the basic theory behind 
Lie-Poisson brackets in Section 2.1, we will show some explicit examples in 
Section 2.2. We then describe general Lie algebra extensions in Section 2.3. 

2.1 Lie-Poisson Brackets on Duals of Lie Algebras 

We begin by taking the Lie algebra q associated with some Lie group. The 
Lie group might be chosen to reflect the symmetries of a physical system. 
There will be a Lie bracket [ , ] : x g ^ g associated with g. Consider the 
dual Q* of 5 with respect to the pairing ( , ) : g* x g — M. Then for real- valued 
functionals F and G, that is, F, G : g* — > M, and ^ e g*, we can define 



The sign choice comes from whether we are considering right invariant (-I-) or 
left invariant (— ) functions on the cotangent bundle of the Lie group [34,9], 
but for our purposes we simply choose the sign as needed. The functional 
derivative is defined by 




(2.1) 



e=0 




(2.2) 
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We shall refer to the bracket [ , ] as the inner bracket and to the bracket { , } 
as the Lie-Poisson bracket. The dual g* together with the Lie-Poisson bracket 
is a Poisson manifold; that is, the bracket { , } is a Lie algebra structure on 
real- valued functionals that is a derivation in each of its arguments. For finite- 
dimensional groups, Eq. (2.1) was first written down by Lie [38] and was 
rediscovered by Berezin [39]; it is also closely related to work of Arnold [40], 
Kirillov [41], Kostant [42], and Souriau [43]. 

The bracket in g is the same as the adjoint action of g on itself: [ o; , /3 ] = ado, P: 
where a, (3 E g. Prom this we define the coadjoint action adj^ of g on g* by 

(adU,/3):=(e, ado/3), (2.3) 

where ^ E g*. We also define the coadjoint bracket [ , ]^ : x g* — > g* by 
[a ,^ ]^ := ad)j, ^, so that 

([a,e]^ ,/3):=(e, [«,/?]); (2.4) 

the bracket [ , ]^ satisfies the identity 

([a,e]^ ,/9) = -([/9,e]^ ,«)■ 



Since the inner bracket is Lie, it satisfies the Jacobi identity, and consequently 
the form given by (2.1) for the Lie-Poisson bracket will automatically satisfy 
the Jacobi identity [44, p. 614]. 

Given a Hamiltonian if : — > IR, the equation of motion for ^ e is 



A 



- e 



, A 



5H 
= T 



- e 



(2.5) 



where A is a Kronecker or Dirac delta, or a combination of both for an infinite- 
dimensional system of several fields. 



2.2 Examples of Lie-Poisson Systems 



We will say that a physical systems can be described by a given Lie-Poisson 
bracket and Hamiltonian if its equations of motion can be written as (2.5); 
the system is then said to be Hamiltonian of the Lie-Poisson type. We give 
four examples: the first is finite-dimensional (the free rigid body. Section 2.2.1) 
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and the second infinite-dimensional (Euler's equation for tlie ideal fluid, Sec- 
tion 2.2.2). The third and fourth examples are also infinite-dimensional and 
serve to introduce the concept of extension. They are low-beta reduced mag- 
netohydro dynamics (MHD) in Section 2.2.3 and compressible reduced MHD 
in Section 2.2.4. These last two examples are meant to illustrate the physical 
relevance of Lie algebra extensions. 



2.2.1 The Free Rigid Body 

The classic example of a Lie-Poisson bracket is obtained by taking for g the 
Lie algebra of the rotation group 5*0(3). If the e(j) denote a basis of = so(3), 
the Lie bracket is given by 

e(i) ,e(j) 

where the c^^ — Sijk are the structure constants of the algebra, in this case 
the totally antisymmetric symbol. Using as a pairing the usual contraction 
between upper and lower indices, with (2.1) we are led to the Lie-Poisson 
bracket 

Sf 1 - k , 9f dg 

where the three- vector £ is in g*, and we have chosen the minus sign in (2.1). 
The coadjoint bracket is obtained using (2.3), 

If we use this coadjoint bracket and insert the Hamiltonian 

in (2.5) we obtain 

Notice how the moment of inertia tensor I plays the role of a metric — it allows 
us to build a quadratic form (the Hamiltonian) from two elements of g*. If 
we take / = diag(/i, I2, 13), we recover Euler's equations for the motion of the 
free rigid body 

and cyclic permutations of 1,2,3. The £i are the angular momenta about the 
axes and the Jj are the principal moments of inertia. This result is naturally 
appealing because we expect the rigid body equations to be invariant under 
the rotation group, hence the choice of SO{3) for G. 
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2.2.2 The Two-dimensional Ideal Fluid 



Consider now an ideal fluid with the flow taking place over a two-dimensional 
domain Jl. Let g be the infinite-dimensional Lie algebra associated with the Lie 
group of volume- preserving diffeomorphisms of In two spatial dimensions 
this is the same as the group of canonical transformations on Q. The bracket 
in is the canonical bracket 

. da dh db da 

dx dy dx dy' 

We formally identify q and Q* and use as the pairing ( , ) the usual integral 
over the fluid domain, 



{F,G)= [ F(x)G(x)d^x, 
Jn 



where x := {x,y). For infinite-dimensional spaces, there are functional ana- 
lytic issues about whether we can make this identification, and take q** = q. 
We will assume here that these relationships hold formally. Sec Marsden and 
Weinstein [32] for references on this subject and Audin [45] for a treatment of 
the identification of q and q*. 

Assuming appropriate boundary conditions for simplicity, we get [ , ]^ = — [ , ] 
from (2.4). (Otherwise the coadjoint bracket would involve extra boundary 
terms.) Take the vorticity lj as the field variable ^ and write for the Hamilto- 
nian 



H[u;]^-Uu,W-'u;), 



where 



( cu) (x) := ^ i<:(x|x') a;(x') d'x', 

and K is Green's function for the Laplacian. The Green's function plays the 
role of a metric since it maps an element of g* (the vorticity uj) into an element 
of Q to be used in the right slot of the pairing. This relationship is only weak: 
the mapping K is not surjective, and thus the metric cannot formally inverted 
(it is called weakly nondegenerate) . When we have identified Q and Q* we shall 
often drop the comma in the pairing and write 

i/[a;] = -|(^0) = i(|V0| 



where uo = defines the streamfunction 0. We work out the evolution 
equation for cu explicitly: 
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u(^) = {uj,H} = [ uj(x.') 
Jn 



5cj(x) SH 



= / a;(x')[5(x-x'),-0(x')] dV 

= y 5(x-x')[cu(x'),0(x')] dV 
= [^(x)>(x)] . 



This is Euler's equation for a two-dimensional ideal fluid. We could also have 
written this result down directly from (2.5) using [ , ]^ = — [ , ]■ 



2.2.3 Low-beta Reduced MHD 

This example will illustrate the concept of a Lie algebra extension, the main 
topic of this paper. Essentially, the idea is to use an algebra of n-tuples, 
which we call an extension, to describe a physical system with more than 
one dynamical variable. As in Section 2.2.2 we consider a flow taking place 
over a two-dimensional domain Q. The Lie algebra Q is again taken to be that 
of volume preserving diffeomorphisms on f2, but now we consider also the 
vector space V of real-valued functions on Vt (an Abelian Lie algebra under 
addition). The semidirect sum of q and V" is a new Lie algebra whose elements 
are two-tuples {a,v) with a bracket deflned by 



(2.7) 



where a and (3 & q, v and w & V. This is a Lie algebra, so we can use the 
prescription of Section 2.1 to build a Lie-Poisson bracket. 



f 




' 5F SG' 




'SF SG' 




'SG SF' 




Jn 




Suj ' Stl) 




Suj ' S%1) 


)\ 



d?x. 



Let UJ — V^(f), where is the electric potential, ip is the magnetic flux, and 
J = V^V' is the current. (We use the same symbol for the electric field as 
for the streamfunction in Section 2.2.2 since they play a similar role.) The 
pairing used is a dot product of the vectors followed by an integral over the 
fiuid domain (again identifying g and g* as in Section 2.2.2). The Hamiltonian 

H[u;;i;] = y^ (|V0r + |V^P) d'x 
with the above bracket leads to the equations of motion 



^ = [-(/;,(/)] . 
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This is a model for low-beta rcdTiccd MHD [12,46,47]. It is obtained by an 
expansion in the inverse aspect ratio e of a tokamak, with e small. This is 
called low beta since the plasma beta (the ratio of plasma pressure to magnetic 
pressure) is of order e^. With a strong toroidal magnetic field, the dynamics 
are then approximately two-dimensional. 

Benjamin [13] used a system with a similar Lie-Poisson structure, but for 
waves in a density-stratified fluid. Semidirect sum structures are ubiquitous 
in advective systems: one variable (in this example, 0) "drags" the others 
along [21]. 



2.2.4 Compressible Reduced MHD 

In general there are other, more general ways to extend Lie algebras besides 
the semidirect sum. The model derived by Hazeltine et al. [17,48] for two- 
dimensional compressible reduced MHD (CRMHD) is an example. This model 
has four fields, and as for the system in Section 2.2.3 it is also obtained from an 
expansion in the inverse aspect ratio of a tokamak. It includes compressibility 
and finite ion Larmor radius effects. The Hamiltonian is 

H[u},v,p,ij] = ^ / 
2 Jn 

where v is the parallel ion velocity, p is the pressure, and f3i is a parameter 
that measures compressibility. The other variables are as in Section 2.2.3. The 
coordinate x points outward from the center of the tokamak in the horizontal 
plane and y is the vertical coordinate. The motion is made two-dimensional 
by the strong toroidal magnetic field. The bracket we will use is 



\V(l)\^ + v^ + 



{p-2i3;xy 
0; 



(2.8) 



- 



UJ 



5F_ 5G 
Scu ' Scu 

SF SG' 
6lj ' 6p 
' SF SG 
Sp ' Sv 



+ v 



SF SG 
Scu ' Sv 



+ 



+ 



+ 



S£ SG 
Sp ' Sou 
SF SG 
Sv ' Sp 



SF SG 
Sv ' Sou 
S£ SG 
Scu ' Sip 



+ 



SF_ SG 
Sip ' Scu 



(2.9) 



Together this bracket and the Hamiltonian (2.8) lead to the equations 



i:u=[u,(P] + [ij,.J] + 2[p,x] 
v^[v,(l)] + [iP,p] + 2(3i[x,ilj] 
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which reduce to the example of Section 2.2.3 in the hmit v—p — Pi — (when 
compressibihty effects are unimportant). 

It is far from clear that the Jacobi identity for (2.9) is satisfied. A direct 
verification is straightforward (if tedious), but we shall see in Section 2.3 that 
there is an easier way. 

2.3 General Algebra Extensions 

We wish to generahze the types of bracket used in Sections 2.2.3 and 2.2.4. 
We build an algebra extension by forming an n-tuple of elements of a single 
Lie algebra g, 

a («!,..., «„), (2.10) 

where «i G g. The most general bracket on this n-tuple space obtained from 
a linear combination of the one in q has components 

n 

[«,/5]a= E ^A^M«/.,/5.] , A=l,...,n, (2.11) 

where the Wx^" are constants. (From now on we will assume that repeated 
indices are summed unless otherwise noted.) Since the bracket in g is anti- 
symmetric the VF's must be symmetric in their upper indices, 

Wx""" = Wx""^ . (2.12) 

This bracket must also satisfy the Jacobi identity 

[/5,7]]a + [/^' [7,a]]A + [7, [«,/3]]a = 0' A = l,...,n. 
The first term can be written 

[a, [/3,7]], = W^A'^^W^.^n«r, [f3,,7.]], 
which when added to the other two gives 

Wx^^W,^" {[ar, [(3,n.]] + [f3r, [7^,«.]] + [7r, [a^,/3.]])=0. 

Wc cannot yet make use of the Jacobi identity in g: the subscripts of a, (3, 
and 7 are different in each term so they represent different elements of g. We 
first relabel the sums and then make use of the Jacobi identity in g to obtain 

{Wx"^ Wr - Wx""" W/^) [ a. , [ , 7. ] ] 

+ {Wx""^ - Wx'"" W/^) , [7. , a. ] ] = . 
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This identity is satisfied if and only if 

Wx"^ W^^"" = Wx"'' W/^ , (2.13) 

which together with (2.12) imphes that the quantity Wx"^^ W„^'^ is symmetric 
in all three free upper indices. If we write the W^s as n matrices W^""^ with 
rows labeled by A and columns by /x, 

[m^M]^" -ly/'^, (2.14) 

then (2.13) says that those matrices pairwise commute: 

^r{y) ^ ^r{'J) (2.15) 

Equations (2.12) and (2.15) form a necessary and sufficient condition: a set 
of n commuting matrices of size n x n satisfying the symmetry given by (2.12) 
can be used to make a good Lie algebra bracket. From this Lie bracket we can 
build a Lie-Poisson bracket using the prescription of (2.1) to obtain 

{F,G}i(0 = ± E wx^'le. 



We now return to the two extension examples of Sections 2.2.3 and 2.2.4 and 
examine them in light of the general extension concept introduced here. 



2. 3. 1 Low- beta Reduced MHD 

For this example we have (^°, ^^) = {u, ■0), with 



(l 


0^ 










.0 






.1 





The reason why we start labehng at will become clearer in Section 4.4. 
The two VF(^) must commute since 1^^°^ = /, the identity. The tensor W 
also satisfies the symmetry property (2.12). Hence, the bracket is a good Lie 
algebra bracket. 



8F SG 
WW 
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2.3.2 Compressihle Reduced MHD 



We have n 

by 



4 and take (^°, {'^) = (a;, f ■?/'), so the tensor W is given 



/ 1 









1 













1 













1 



-AO y 





1 

















1^(3) 



-A y 

^ 




1 y 



(2.16) 



It is easy to verify that these matrices commute and that the tensor W satisfies 
the symmetry property, so that the Lie-Poisson bracket given by (2.9) satisfies 
the Jacobi identity. (See Section 4.4 for an explanation of why the labeling is 
chosen to begin at zero.) 



3 Extension of a Lie Algebra 



In this section we review the theory of Lie algebra cohomology and its ap- 
plication to extensions. This is useful for shedding light on the methods used 
in Section 4 for classifying extensions. However, the mathematical details pre- 
sented in this section can be skipped without seriously compromising the fiavor 
of the classification scheme of Section 4. 



3.1 Cohomology of Lie Algebras 



We now introduce the abstract formalism of Lie algebra cohomology. Histori- 
cally there were two different reasons for the development of this theory. One, 
known as the Chevalley-Eilenberg formulation [27], was developed from de 
Rham cohomology. de Rham cohomology concerns the relationship between 
exact and closed differential forms, which is determined by the global prop- 
erties (topology) of a differentiable manifold. A Lie group is a differentiable 
manifold and so has an associated de Rham cohomology. If invariant differen- 
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tial forms are used in the computation, one is led to the cohomology of Lie 
algebras presented in this section [28,29,49]. The second motivation is the one 
that concerns us: we will show in Section 3.2 that the extension problem — the 
problem of enumerating extensions of a Lie algebra — can be related to the 
cohomology of Lie algebras. 

Let be a Lie algebra, and let the vector space V over the field K (which we 
take to be the real numbers later) be a left g-module,[^ that is, there is an 
operator p : g x V ^ V such that 



pa (V + v')=PaV + Pa v' , 
Pa+a' V = PaV + Pa' V, 
P[a,a']V= [Pa,Pa'] V , (3.1) 

for a, a' G g and v, v' G V . The operator p is known as a left action. A g- 
module gives a representation of g on V . 

An n-dimensional l^-valued cochain ujn for g, or just ra-cochain for short, is a 
skew-symmetric n-linear mapping 

c^„:gxgx---xg — >V. 

Cochains are Lie algebra cohomology analogues of differential forms on a man- 
ifold. Addition and scalar multiplication of n-cochains are defined in the ob- 
vious manner by 



(cu„ + u;^)(ai, . . . , a„) :=cj„(ai, ...,«„)+ • • • , "n), 

(au;„)(ai, . . . , a„) :=at^„(ai, . . .,«„), 



where ai, . . . , a„ G g and a G K. The set of all n-cochains thus forms a vector 
space over the field K and is denoted by C"(g, V). The 0-cochains are defined 
to be just elements of V, so that C°(g, V) = V. 

The coboundary operator is the map between cochains. 



5n:C"(g,V)— C"+i(g,r), 



defined by 



^ When y is a right g-module, we have P[a,a'] = ~[pa ,Pa']- The results of this 
section can be adapted to a right action by changing the sign every time a commu- 
tator appears. 
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n+1 

{SnU;n){ai, . . . ,an+l) := Xl(-)'^Vai'^n(ttl, • • • , Qfi, • • • , "n+l) 

i=l 

n+1 

+ i-y^''^n{[otj ,oik] ,ai, . . . ,aj, . . . ,ak, . . • ,Q;n+i), 

j,k=l 
3<k 

where the caret means an argument is omitted. We shall often drop the n 
subscript on s„, deducing it from the dimension of the cochain on which s 
acts. 

We shall make use mostly of the first few cases 

(scJo)(ai) = Pai ^0, (3.2) 

(sc^l)(ai,a2) =Pai t^l(tt2) - Pa2 t^l(«l) " ^^1 ( [ "l , ^2 ] ) , (3.3) 
(su;2)(ai, 02, tta) =Pai '^2(a2, "3) + Pa2 '^2(a3, Oil) + Pag '^2(ai, "2) 

-a;2([Q;i ,0:2] ,0:3) - a;2([Q;2 , as] , ^i) - a;2([a3 , ] , 0:2) ■ (3.4) 

It is easy to verify that sojn defines an (n + l)-cochain, and it is straight- 
forward (if tedious) to show that s„+is„ ~ — Q. For this to be true, The 
homomorphism property of p is crucial. 

An n-cocycle is an element Un of C'^{q,V) such that SnOJn — 0. An n- 
coboundary cjcob is an element of C"(g, K) for which there exists an ele- 
ment ujn-i of C"~^(g, V) such that a;cob = sujn-i- Note that all coboundaries 
are cocycles, but not vice-versa. 

Let 

Z;(0,y) :=kers„ 

be the vector subspace of all n-cocycles, (3, V) C C"(0, V), and let 

-Sp(0,l^) — range Sn-i 

be the vector subspace of all n-coboundaries, B'^{q^V) C C"(g, K). The nth 
cohomology group of q with coefficients in V is defined to be the quotient 
vector space 

h;{q,v):^z;{q,v)/b;{q,v). (3.5) 

Note that for n > dimg, we have i/^(0, V) = Z^(0, V) = ^^(g, V) = 0. 
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3.2 Application of Cohomology to Extensions 

In Section 2.3 we gave a definition of extension that is specific to our problem. 
We will now define extensions in a more abstract manner. We then show 
how the cohomology of Lie algebras of Section 3.1 is related to the problem of 
classifying extensions. In Section 4 we will return to the more concrete concept 
of extension, of the form given in Section 2.3. 

Let /i : 0i — > Qi+i be a collection of Lie algebra homomorphisms, 

i *- 01+1 *- 01+2 s- ■ ■ ■ . 

The sequence fi is called an exact sequence of Lie algebra homomorphisms if 

range fi = ker /j+i . 

Let Q, f), and o be Lie algebras. The algebra [) is said to be an extension of g 
by a if there is a short exact sequence of Lie algebra homomorphisms 

-a^-{)^^0 -0. (3.6) 

The homomorphism i is an insertion (injection), and tt is a projection (sur- 
jection). We shall distinguish brackets in the different algebras by appropriate 
subscripts. Wc also define r : g ^ [) to be a linear mapping such that tt o r = l|g 
(the identity mapping in q). Note that r is not unique, since the kernel of n 
is not trivial. Let P e i), rj E a; then 

using the homomorphism property of tt and ir oi — 0, a consequence of the 
exactness of the sequence. Thus [P ,ir)]f^ e ker tt = range i, and i a is an ideal 
in f) since ,irj] e ia. Hence, we can form the quotient algebra {)/a, with 
equivalence classes denoted by P + a. By exactness 7r(/3 + a) = tt/S, so g is 
isomorphic to l)/a and we write g = [)/a. 

Though i a is a subalgebra of f), rg is not necessarily a subalgebra of 1^, for in 
general 

[ra,T/5](, ^ r 

for G 0; that is, r is not necessarily a homomorphism. The classification 
problem essentially resides in the determination of how much r differs from 
a homomorphism. The cohomology machinery of Section 3.1 is the key to 
quantifying this difference, and we proceed to show this. 
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To this end, we use the algebra a as the vector space V of Section 3.1, so 
that o will be a left g-module. We define the left action as 



(3.7) 



for a E Q and r/ G a. For a to be a left g-module, we need p to be a homomor- 
phism, i.e., p must satisfy (3.1). Therefore consider 



^Pai'^[T(3,ir)]f^- ppi~^[ra,ir]]f^ 



to;, [T/3,i77]( 



— I 



TP,[Ta,iri]f 



which upon using the Jacobi identity in [) becomes 



([ra,r/3]f, -r ,iri 



(3.8) 



By applying n on the expression in parentheses of the last term of (3.8), we 
see that it vanishes and so is in ker tt, and by exactness it is also in i a. Thus 
the f) commutator above involves two elements of i a. We define a; : g x — > a 

by 

u{a,P) :=t-'{[Ta,rf3]^-T[a,(3]^). (3.9) 
The mapping is well defined on i a. Equation (3.8) becomes 

[Pa,Pp]V = Pla,p]^V+[^{<^,P) ,V]a- (3-10) 

Therefore, p satisfies the homomorphism property if either of the following is 
true: 

(i) a is Abehan, 

(ii) T is a homomorphism. 

Condition (i) implies [ , = 0, while condition (ii) means 

[tq;,t/3][, ^T[a,(3]g, 

which implies = 0. If either of these conditions is satisfied, a with the 
action p is a left g-module. We treat these two cases separately in Sections 3.3 
and 3.4, respectively. 
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3.3 Extension by an Abelian Lie Algebra 



In this section we assume that the homomorphism condition (i) at the end of 
Section 3.2 is met. Therefore o is a left g-module, and we can define a-valued 
cochains on g. In particular, uj defined by (3.9) is a 2-cochain, uj G C^(g, a), 
that measures the "failure" of r to be a homomorphism. We now show, more- 
over, that a; is a 2-cocycle, uj G Z'^{g, a). By using (3.4), 



-'^([«>/5]0>7) -^([/5,7]0,«) -'^([7,«]g,/?) , 



ra, [t/3,t7] 



+ eye. perm. 



+ eye. perm. 



0. 



The first parenthesis vanishes by the Jacobi identity in 1^, the second by the 
Jacobi identity in g, and the other terms were canceled in pairs. Hence, cu is 
a 2-cocycle. 

Two extensions 1^ and f)' are equivalent if there exists a Lie algebra isomor- 
phism (7 such that the diagram 



a 







(3.11) 



is commutative, that is if cr o i = i' and t: — tt' o a. 

There will be an injection r associated with tt and a r' associated with tt', such 
that TT o r = l|g = vr' o r'. The linear map u = a~^T' — t must be from g to i a, 
so G C^(g, o). Consider p and p' respectively defined using r, i and T',i' 
by (3.7). Then 



^i~^[T a ,ir]\- i-^[{v + t) a ,ir} \ , 
= -r^[i/Q;,i77]^ = 0, 



(3.12) 



since a is Abelian. Hence r and r' define the same p. Now consider the 2- 
cocycles u and u' defined from r and r' by (3.9). We have 
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= Pa P) - Pn «) - i'^i^ [ a , ]g- 

Comparing this with (3.3), we see that 

uj' -uj^s{i-^v), (3.13) 

so uj and uo' differ by a coboundary. Hence they represent the same element 
in Hp{Q,a). Equivalent extensions uniquely define an element of the second 
cohomology group H'^{q, a). Note that this is true in particular for f) = ()', cr = 
1, so that the element of Hp{g, o) is independent of the choice of r. 

We are now ready to write down explicitly the bracket in f). We can repre- 
sent an element a e 1) as a two-tuple: a = (cti, where ai G g and 0:2 G a 
(J) = © ci as a vector space). The injection i is then ia2 — (0, 02), the projec- 
tion TT is TT (cKi, 0:2) = «!, and since the extension is independent of the choice 
of T we take rai — (cvi, 0). By linearity 

[a, = [ (ai, 0), (A, 0) ]^ + [ (0, a2), (0, /^a) ]^ 

+ [(ai,0),(0,/32)](, + [(0,a2),(/3i,0)](,. 

We know that [(0,q;2), (0, /32)][, = since is Abehan. By definition of the 
cocycle a;, Eq. (3.9), we have 

[(ai,0),(/3i,0)](,= [rai,rA]^ 

Finally, by the definition of p, Eq. (3.7), 

[(ai,0), (0,/32)](, = [Tai,i/32](, = Pa^ P2, 
and similarly for [ (0, 0:2), (A) 0) ][^, with opposite sign. So the bracket is 

[«>/5]f,= Pai/52-P/3ia2+a;(ai,/3i)^. (3.14) 

As a check we work out the Jacobi identity in f): 
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« J/5>7][,]^ = ([«! , [/^^Tlilg, Pai [/5>7]2 -P[/3,7]i Q;2 +i^(a;i, [/5,7]i)) 

= (["l'[A'7l]fl]g, Pai(P/3i72 -P7i/52+^(A,7l)) 

- ,71 ], «2 + [A , 7l ]g)) ■ 

Upon adding permutations, the first component will vanish by the Jacobi 
identity in q. We are left with 



a 



J/5'7](,](j + CyC. perm.= (^0, (paiP/3i - Pp^Pa^ - P[ai 72 

+ Pai^^(/5i,7i) -^(["i,A]0,7i)) +cyc. perm.. 



which vanishes by the the homomorphism property of p and the fact that uj 
is a 2-cocycle, Eq. (3.4). 

Equation (3.14) is the most general form of the Lie bracket for extension by 
an Abelian Lie algebra. It turns out that the theory of extension by a non- 
Abehan algebra can be reduced to the study of extension by the center of a, 
which is Abehan [29]. We will not need this fact here, as the only extensions by 
non- Abelian algebras we will deal with are of the simpler type of Section 3.4. 

We have thus shown that equivalent extensions are enumerated by the second 
cohomology group H'^{q, a). The coordinate transformation a used in (3.11) to 
define equivalence of extensions preserves the form of g and a as subsets of [). 
However, we have the freedom to choose coordinate transformations which do 
transform these subsets. All we require is that the isomorphism a between 1^ 
and i)' be a Lie algebra isomorphism. We can represent this by the diagram 

>a^-f)^^0 -0 (3.15) 



a' ^ f)' ^ S' 0. 

The primed and the unprimed extensions are not equivalent, but they are 
isomorphic [50] . Cohomology for us is not the whole story, since we are inter- 
ested in isomorphic extensions, but it will guide our classification scheme. We 
discuss this point further in Section 4.3. 



3.4 Semidirect and Direct Extensions 



Assume now that uj defined by (3.9) is a coboundary. By (3.13) there exists 
an equivalent extension with cu = 0. For that equivalent extension, r is a ho- 
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momorphism and condition (ii) at the end of Section 3.2 is satisfied. Thus the 
sequence 



^-^0^ 

is an exact sequence of Lie algebra homomorphisms, as well as the sequence 
given by (3.6). We then say that the extension is a semidirect extension (or a 
semidirect sum of algebras) by analogy with the group case. More generally, 
we say that f) splits if it is isomorphic to a semidirect sum, which corresponds 
to uj being a coboundary, not necessarily zero. If a is not Abelian, then (3.12) 
is not satisfied and two equivalent extensions (or two different choices of r) 
do not necessarily lead to the same p. 

Representing elements of f) as 2-tuples, as in Section 3.3, we can derive the 
bracket in f) for a semidirect sum, 

([«l>/^l]fl> Pai/32-P/3ia2 + [a2,/?2]a), (3-16) 

where we have not assumed a Abelian. Verifying Jacobi for (3.16) we find 
the p must also satisfy 

which is trivially satisfied if a is Abelian, but in general this condition states 
that Pa is a derivation on o. 

Now consider the case where i^^ is a homomorphism and keri~^ = ranger. 
Then the sequence 

is exact in both directions and, hence, both i and tt = are bijections. The 
action of on a is 







since by exactness i^^ o r = 0. This is called a direct sum. Note that in this 
case the role of g and a is interchangeable and they are both ideals in [). The 
bracket in f) is easily obtained from (3.16) by letting p — 0, 

[a,/3]^=([ai,/3i]g,[a2,/32]„). (3.17) 

Semidirect and direct extensions play an important role in physics. A simple 
example of a semidirect extension structure is when q is the Lie algebra so(3) 
associated with the rotation group 5*0(3) and a is M^. Their semidirect sum 



19 



is the algebra of the six parameter EucUdean group of rotations and trans- 
lations. That algebra can be used in a Lie-Poisson bracket to describe the 
dynamics of the heavy top (see for example [11,51]). We have already dis- 
cussed the semidirect sum in Section 2.2.3. The bracket (2.7) is a semidirect 
sum, with Q the algebra of the group of volume-preserving diffeomorphisms 
and a the Abelian Lie algebra of functions on M^. The action is just the adjoint 
action paV :— [a ,v] obtained by identifying g and a. 

A Lie-Poisson bracket built from a direct extension is just a sum of the sep- 
arate brackets. The interaction between the variables can only come from 
the Hamiltonian or from constitutive equations. For example in the baro- 
clinic instability model of two superimposed two fluid layers with different 
potential vorticities the two layers are coupled through the potential vorticity 
relation [14]. 



4 Classification of Extensions of a Lie Algebra 



In this section we return to the main problem introduced in Section 2.3: 
the classification of algebra extensions built by forming n-tuples of elements 
of a single Lie algebra q. The elements of this Lie algebra f) are written 
as a := {ai, . . . , «„), ctj e q, with a bracket defined by 

[«,/?], = W^a'^'^ [a,,P,], (2.11) 

where W^'^" are constants. We will call n the order of the extension. Recall 
(see Section 2.3) the VF's are symmetric in their upper indices, 

Wx"" ^ Wx"" , (2.15) 

and commute, 

lyH^yW = (2.15) 

where the n x n matrices W'^"^ are defined by [W'^''^^'' := Wx"^. Since the W^s 
are 3-tensors we can also represent their elements by matrices obtained by 
fixing the lower index. 



W, 



w, 



(A) 



Wx''", (4.1) 



which are symmetric but do not commute. Either collection of matrices, (2.14) 
or (4.1), completely describes the Lie bracket, and which one we use will be 
understood by whether the parenthesized index is up or down. 

What do we mean by a classification? A classification is achieved if we obtain 
a set of normal forms for the extensions that are independent, that is, not 
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related by linear transformations. We use linear transformations because they 
preserve the Lie-Poisson structure — they amount to transformations of the W 
tensor. We thus begin by assuming the most general W possible. 

We first show in Section 4.1 how an extension can be broken down into a direct 
sum of degenerate subblocks (degenerate in the sense that the eigenvalues have 
multiplicity greater than unity). The classification scheme is thus reduced to 
the study of a single degenerate subblock. In Section 4.2 we couch our par- 
ticular extension problem in terms of the Lie algebra cohomology language of 
Section 3.2 and apply the techniques therein. The limitations of this cohomol- 
ogy approach are investigated in Section 4.3, and we look at other coordinate 
transformations that do not necessarily preserve the extension structure of the 
algebra, as expressed in diagram (3.15). In Section 4.5 we introduce a partic- 
ular type of extension, called the Leibniz extension, which is in a sense the 
"maximal" extension. Finally, in Section 4.6 we give an explicit classification 
of solvable extensions up to order four. 

4-1 Direct Sum Structure 

A set of commuting matrices can be put into simultaneous block-diagonal form 
by a coordinate transformation, each block corresponding to a degenerate 
eigenvalue [52]. Let us denote the change of basis by a matrix M^", with 
inverse (M~^)-^, such that the matrix VF^^^ whose components are given by 

is in block-diagonal form for all i/ [52] . However, Wx^" is a 3-tensor and so the 
third index is also subject to the coordinate change: 

This last step adds linear combinations of the W^'^'^^s together, so the W^'^'^^s 
and the W^^^^s have the same block-diagonal structure. Note that the W ten- 
sors are still symmetric in their upper indices, since this property is preserved 
by a change of basis. So from now on we just assume that we are working 
in a basis where the W^'^^^s are block-diagonal and symmetric in their upper 
indices; this symmetry means that if we look at a 1^ as a cube, then in the 
block-diagonal basis it consists of smaller cubes along the main diagonal. This 
is the 3-tensor equivalent of a block-diagonal matrix. 

Block-diagonalization is the first step in the classification: each block of W is 
associated with an ideal (hence, a subalgebra) in the full n-tuple algebra q. 
Hence, by the definition of Section 3.4 the algebra g is a direct sum of the sub- 
algebras associated with each block. Each of these subalgebras can be studied 
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independently, so from now on we assume that we have n commuting matri- 
ces, each with n-fold degenerate eigenvalues. The eigenvalues can, however, be 
different for each matrix. 



Such a set of commuting matrices can be put into lower-triangular form by a 
coordinate change, and again the transformation of the third index preserves 
this structure (though it changes the eigenvalue of each matrix). The eigen- 
value of each matrix lies on the diagonal; we denote the eigenvalue of W^^^ 
by A^'*). The matrix VF^fi), which as prescribed by (4.1) consists of the first row 
of the lower-triangular matrices W'^^\ is given by 



^ A« ••• 0^ 
A(2) • • • 



^A(") ••• Oj 



Evidently, the symmetry of requires 

AM 

that is, all the matrices W^^'^ are nilpotent (their eigenvalues vanish) except 
for VT'^^^ when 6' 7^ 0. If this first eigenvalue is nonzero then it can be scaled 
to ^ = 1 by the coordinate transformation Mj," = 9'^ Si,"'. We will use the 
symbol 9 to mean a variable that can take the value or 1. 



4-2 Connection to Cohomology 



We now bring together the abstract notions of Section 3 with the n-tuple 
extensions of Section 2.3. It is shown in Section 4.2.1 that we need only classify 
the case of ^ = 0. This case will be seen to correspond to solvable extensions, 
which we classify in Section 4.2.2. 



4.2.1 Preliminary Splitting 

Assume wc arc in the basis described at the end of Section 4.1 and, for now, 
suppose 9 = 1. The set of elements of the form P = {0, P2, ■ ■ ■ , Pn) is a nilpotent 
ideal in f) that we denote by a (a is thus a solvable subalgebra [53]). Hence, we 
can construct the algebra Q = f)/a, so that is an extension of q by a. If is 
semisimple, then is the radical of f) (the maximal solvable ideal). It is easy 
to see that the elements of Q are of the form a = (ai, 0, . . . , 0). We will now 
see that [) splits; that is, there exist coordinates in which [) is manifestly the 
semidirect sum of g and the (in general non-Abelian) algebra a. 
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In Appendix A we give a lower-triangular coordinate transformation that 
makes W^^^ = I, the identity matrix. Assuming we have effected this trans- 
formation, the mappings i, vr, and r of Section 3.2 are given by 

i : a — > [), i{a2, ...,«„) = (0, 0^2, ... , an) 
TT : [) — > g, 7r(ai, a2, . . . ,an) = ai, 
t:q — >i), r(ai) = (ai,0,...,0), 

and the cocycle of Eq. (3.9) is 

= [(«i,0,...,0),(A,0,...,0)]^-([«i,A],0,...,0) 
= 0. 

Since uj = 0, the extension is a semidirect sum (see Section 3.4). The coor- 
dinate transformation that made W^^^ = I removed a coboundary, making 
the above cocycle vanish identically. For the case where g is finite- dimensional 
and semisimple, we have an explicit demonstration of the Levi decomposition 
theorem: any finite-dimensional|^ Lie algebra f) (of characteristic zero) with 
radical a is the semidirect sum of a semisimple Lie algebra g and a [53] . 



4-2.2 Solvable Extensions 

Above we assumed the eigenvalue 6 of the first matrix was unity; however, if 
this eigenvalue vanishes, then we have a solvable algebra of n-tuples to begin 
with. Since n is arbitrary we can study these two solvable cases together. 

Thus, we now suppose f) is a solvable Lie algebra of n-tuples (we reuse the 
symbols f), g, and a to parallel the notation of Section 3.1), where all of the 
theiy^'^^'s are lower-triangular with zeros along the diagonal. Note that W^"^^ = 
0, so the set of elements of the form a = (0, . . . , 0, an) forms an Abelian sub- 
algebra of f). In fact, this subalgebra is an ideal. Now assume f) contains an 
Abelian ideal of order n — m (the order of this ideal is at least 1), which we 
denote by a. The elements of a can always be cast in the form 

a= (0, . . . , 0, ftm+i, •••,«„) 

via a coordinate transformation that preserves the lower-triangular, nilpotent 
form of the W^^^'' . 



^ The inner bracket can be infinite dimensional, but the order of the extension is 
finite. 
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We also denote by Q the algebra of m-tuples with the bracket 

m 

[(«!,..., «„),(/?!,..., ^ Wx^" [a,, , (3^] , X=l,...,m. 

IJ,,u=l 

It is trivial to show that g = [)/a, so that i) is an extension of q by a. Since o is 
Abelian we can use the formalism of Section 3.1 (the other case we used above 
was for non- Abelian but where the extension was semidirect). The injection 
and projection maps are given by 

i : a — > [), . . . , a„) = (0, . . . , 0, ttm+i, ■ ■ • , "n), 

T : — > [), r(ai, .... am) = («!, • • • , a™, 0, . . . , 0). 
Prom the definition of the action, Eq. (3.7), we have for a e g and r] & a, 



= [ (ai, . . . , 0, . . . , 0) , (0, . . . , 0, rj^+i, ■ ■ ■ , ?7n) ](, 

m n— 1 

= E E iO,...,0,Wj:i2[c^^,Vu],...,W^'"'[a,,v.]). (4.2) 

IJ,=1 v=m+l 

In addition to the action, the solvable extension is also characterized by the 
cocycle defined in Eq. (3.9), 



= [ («!, . . . , a^, 0, . . . , 0) , (/3i, . . . , 0, . . . , 0) ](, 

- r [ («!, . . . , am) , iPl, ■■■,Pm)]g 

m 

= E (0,...,0,WC+l[«/.,/5.],---,W^n'^l«/.,/3.]). 



(4.3) 



We can illustrate which parts of the W's contribute to the action and which 
to the cocycle by writing 



Wa 


ta 


T 






A = m + 1, . . . , n. 



(4.4) 



where the wa's are m x m symmetric matrices that determine the cocycle uu 
and the Ta's are m x {n — m) matrices that determine the action p. The 
{n — m) X {n — m) zero matrix on the bottom right of the VF(a)'s is a con- 
sequence of being Abehan. 

The algebra g is completely characterized by the I^(a), A = 1, . . . ,m. Hence 
we can look for the maximal Abelian ideal of g and repeat the procedure we 
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used for the full t). It is straightforward to show that although coordinate 
transformations of g might change the cocycle cu and the action p, they will 
not alter the form of (4.4) . 

Recall that in Section 3.1 we defined 2-coboundaries as 2-cocycles obtained 
from 1-cochains by the coboundary operator, s. The 2-coboundaries turned 
out to be removable obstructions to a semidirect sum structure. Here the 
coboundaries are associated with the parts of the W(^x) that can be removed 
by (a restricted class of) coordinate transformations, as shown below. 

Let us explore the connection between 1-cochains and coboundaries in the 
present context. Since a 1-cochain is just a linear mapping from q to a, 
for a — (cti, . . . , am) G g we can write this as 

m 

•^^H^) = - ^/ , /x = m-M,...,n, (4.5) 

where the k^^ arc arbitrary constants. To find the form of a 2-coboundary we 
act on the 1-cochain (4.5) with the coboundary operator; using (3.3) and (4.2) 
we obtain 



= p«^W(/3)+p^^W(a)-^W([«,/?]J, 



m n 



/i=l v=m+\ 



(4.6) 



m n 



//=1 v=m+\ ix,v,cr=l 

After inserting (4.5) into (4.6) and relabeling, we obtain the general form of 
a 2-coboundary 



where 



^°b(a,/?)= ^A^n«/.,/5.], A = m+l,...,n, 

IJ,,v=\ 



m n 

^'''^^^E^a'W^/'^- E (ka^'Wr + K-'Wr). (4.7) 

T=l a=m+l 



To see how coboundaries are removed, consider the lower-triangular coordinate 
transformation 








k 


cl 
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where a labels rows. This transformation subtracts V(a) from ^^(a) for X > m 
and leaves the first m of the H^(a)'s unchanged. In other words, if W is the 
transformed W, 



w. 



(A) , 



(WA - Va) 



X — 1, ... ,171; 
, X — m + 1, . . . ,n. 







We have also included in this transformation an arbitrary scale factor c. Since 
by (4.3) the block in the upper-left characterizes the cocycle, we see that the 
transformed cocycle is the cocycle characterized by wa minus the coboundary 
characterized by Va. 

The special case we will encounter most often is when the maximal Abelian 
ideal of f) simply consists of elements of the form (0, . . . , 0, «„). For this 
case m = n — 1, and the action vanishes since Wn^^ = (the extension is 
central). The cocycle uu is entirely determined by W(^n)- The form of the co- 
boundary is reduced to 



n-l 



(4.9) 



T=l 



that is, a linear combinations of the first (n — 1) matrices. Thus it is easy to 
see at a glance which parts of the cocycle characterized M^(n) can be removed 
by lower-triangular coordinate transformations. 



4-3 Further Coordinate Transformations 



In the previous section we restricted ourselves to lower-triangular coordinate 
transformations, which in general preserve the lower-triangular structure of 
the W'^^\ But when the matrices are relatively sparse, there exist non- lower- 
triangular coordinate transformations that nonetheless preserve the lower- 
triangular structure. As alluded to in Section 3.3, these transformations are 
outside the scope of cohomology theory, which is restricted to transformations 
that preserve the exact form of the action and the algebras g and a, as shown 
by (4.8). In other words, cohomology theory classifies extensions given g, a, 
and p. We need not obey this restriction. We can allow non-lower-triangular 
coordinate transformations as long as they preserve the lower-triangular struc- 
ture of the W^^'^'s. 

We now discuss a particular class of such transformations that will be useful 
in Section 4.6. Consider the case where both the algebra of (n — l)-tuples q 
and that of 1-tuples a are Abehan. Then the possible (solvable) extensions, in 
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lower triangular form, are characterized by W^ 



(A) 



0, A = 1 n — 1, with 



W(^n) arbitrary (except for Wn^"' = 0). Let us apply a coordinate change of the 
form 



M 



where m is an (n — 1) x [n — 1) nonsingular matrix and c is again a nonzero 
scale factor. Denoting by W the transformed W, we have 




W 



(A) 















A = l,. 

X = n. 



.n-h 



(4.10) 



This transformation does not change the lower-triangular form of the exten- 
sion, even if m is not lower-triangular. The manner in which w„ is transformed 
by M is very similar to that of a (possibly singular) metric tensor: it can be 
diagonalized and rescaled such that all its eigenvalues are or ±1. We can 
also change the overall sign of the eigenvalues using c (something that cannot 
be done for a metric tensor). Hence, we shall order the eigenvalues such that 
the -fl's come first, followed by the — I's, and finally by the O's. We will show 
in Section 4.6 how the negative eigenvalues can be eliminated to harmonize 
the notation. 



4-4 Appending a Semisimple Part 



In Section 4.2 we showed that because of the Levi decomposition theorem we 
only needed to classify the solvable part of the extension for a given degen- 
erate block. Most physical applications have a semisimple part (^ = 1); when 
this is so, we shall label the matrices by W^^\W^^\ . . . ,W^'^\ where they 
are now of size n + I and W^'^^ is the identity.]^ Thus the matrices labeled 
by W^^\ . . . , W^"^^ will always form a solvable subalgebra. This explains the 
labeling in Sections 2.3.1 and 2.3.2. 

If the extension has a semisimple part (9 = 1, or equivalently W^^^ = /), we 
shall refer to it as semidirect. This was the case treated in Section 4.2.1. If 
the extension is not semidirect, then it is solvable (and contains n matrices 
instead of n + 1). 

^ The term semisimple is not quite precise: if the base algebra is not semisimple then 
neither is the extension. However we will use the term to distinguish the different 
cases. 
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Given a solvable algebra of n-tuples we can carry out in some sense the in- 
verse of the Levi decomposition and append a semisimple part to the ex- 
tension. Effectively, this means that the n x n matrices W^^\ . . . , W^"^ are 
made n-|-lxn-|-lby adding a row and column of zeros. Then we simply 
append the matrix l^^*^) = / to the extension. In this manner we construct 
a semisimple extension from a solvable one. This is useful since we will be 
classifying solvable extensions, and afterwards we will want to recover their 
semi direct counterpart. 

The extension obtained by appending a semisimple part to the completely 
Abelian algebra of n-tuples will be called pure semidirect. It is characterized 
by 1^(0) _ y^^f^u = for /X, 1/ > 0. 



4-5 Leibniz Extension 



A particular extension that we shall consider is called the Leibniz exten- 
sion [54] . For the solvable case this extension has the form 



/ 





1 
1 



(4.11) 



V 1 V 

or Wx'^^ — Sx-i'^, A > 1; i.e. the first matrix is an n x n Jordan block. In this 
case the other matrices, in order to commute with W^^\ must be in striped 
lower-triangular form [52] . After using the symmetry of the upper indices the 
matrices can be reduced to 

^r{^^) ^ (^jyy^ (412) 

where on the right-hand side the u denotes an exponent, not a superscript. 
An equivalent way of characterizing the Leibniz extension is 

W^,u^Sx'+\ i,,u,X^l,...,n. (4.13) 

The tensor 5 is an ordinary Kronecker delta. Note that neither (4.12) nor (4.13) 
are covariant expressions, reflecting the coordinate-dependent nature of the 
Leibniz extension. 

The Leibniz extension is in some sense a "maximal" extension: it is the only 
extension that has W(^x) ior all X — 2, . . . ,n (up to coordinate transfor- 
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mations). Its uniqueness will become clear in Section 4.6, and is discussed in 
Thiffeault [55]. 



To construct the semidirect Leibniz extension, we append W^^^ — I, & square 
matrix of size n + 1, to the solvable Leibniz extension above, as described in 
Section 4.4. 



4-6 Low- order Extensions 

We now classify the algebra extensions of low order. As demonstrated in Sec- 
tion 4.2 we only need to classify solvable algebras, which means that = 
for all cases. We will do the classification up to order n = 4. For each case we 
first write down the most general set of lower-triangular matrices W'^''^ (we 
have already used the fact that a set of commuting matrices can be lower- 
triangularized) with the symmetry W\'^'^ = W\^^ built in. Then we look at 
what sort of restrictions the commutativity of the matrices places on the el- 
ements. Finally, we eliminate coboundaries for each case by the methods of 
Sections 4.2 and 4.3. This requires coordinate transformations, but we usually 
will not bother using new symbols and just assume the transformation was 
effected. 

Note that, due to the lower-triangular structure of the extensions, the classi- 
fication found for an m-tuple algebra applies to the first m elements of an n- 
tuplc algebra, n > m. Thus, W(^n) is the cocycle that contains all of the new 
information not included in the previous m — n — 1 classification. These com- 
ments will become clearer as we proceed. 

We shall call an order n extension trivial if W(^n) = 0, so that the cocycle 
appended to the order n — 1 extension contributes nothing to the bracket. 



4.6.1 n=l 

This case is Abelian, with the only possible element Wi'^ — 0. 



4.6.2 n=2 

The most general lower-triangular form for the matrices is 












29 



If VTa^^ 7^ 0, then we can rescale it to unity. Hence we let W2^^ '■= Oi, 
where = or 1. The case = is the Abehan case, while for 9 = 1 
we have the n = 2 Leibniz extension (Section 4.5). Thus for n = 2 there are 
only two possible algebras. The cocycle which we have added at this stage is 
characterized by 9i. 



4.6.3 n=3 

Using the result of Section 4.6.2, the most general lower-triangular form is 



^ ^ 

^1 
[Ws^'Ws^' ) 



^ ^ 


yw^^'Ws^^ J 



and W^^^ = 0. These satisfy the symmetry condition (2.12), and the require- 
ment that the matrices commute leads to the condition 



0^ ^ 0. 



The symmetric matrix representing the cocycle is 



^1^3 11 1^321 ^ 



W^(3) 



V 



w/321 ^^322 




(4.14) 



If 01 = 1, then W-i'^^ must vanish. Then, by (4.9) we can remove from W(^s) 
a multiple of W(2), and therefore we may assume W^^^ vanishes. A suitable 
rescaling allows us to write W^^^ = 6*2, where 6'2 = or 1. The cocycle for the 
case ^1 = 1 is thus 



^ 62 ^ 



W^(3) 



^2 





For ^2 = 1 we have the Leibniz extension (Section 4.5). 

If 9i = 0, we have the case discussed in Section 4.3. For this case we can 
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diagonalize and rescale 1^(3) such that 



W^(3) = 



Ai 
A2 




where (Ai, A2) can be (1, 1), (1, 0), (0, 0), or (1, —1). This last case, as alluded to 
at the end of Section 4.3, can be transformed so that it corresponds to 61 = 0, 
6*2 = 1. The choice (1,0) can be transformed to the ^1 = 1, 6*2 = case. 
Finally for (Ai, A2) = (1, 1) we can use the complex transformation 



to transform to the ^1 = 0, ^2 = 1 case. 



We allow complex transformations in our classification because we are chiefiy 
interested in finding Casimir invariants for Lie-Poisson brackets. If we dis- 
allowed complex transformations, the final classification would contain a few 
more members. The use of complex transformations will be noted as we pro- 
ceed. 

There are thus four independent extensions for n = 3, corresponding to 

(^i,^2)e{(o,o),(o,i),(i,o),(i,i)}. 

These will be referred to as Cases 1-4, respectively. Cases 1 and 3 have 62 — 0, 
and so are trivial (11^(3) = 0). Case 2 is the solvable part of the compressible 
reduced MHD bracket (Section 2.3.2). Case 4 is the solvable Leibniz extension. 



4.6.4 n—4 

Proceeding as before and using the result of Sections 4.6.2 and 4.6.3, we now 
know that we need only write 



W, 



(4) 



^ W^^ w^^ w/^ 

W^^ W^'^ Wi''^ 

w^^ w/^ w/^ 



V 







'4 









(4.15) 



The matrices ^^(i), W(2), and 14^(3) are given by their n = 3 analogues padded 
with an extra row and column of zeros (owing to the lower-triangular form 
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of the matrices). The requirement that the matrices W^^^ . . . VF^^^ commute 
leads to the conditions 



02 W4^^ 


= 0, 




= eiW4 




= 0, 


01 W4^^ 


= 0. 



22 



(4.16) 



There are four cases to look at, corresponding to the possible values of 9i 
and 6*2. 

Case 1 01^0,02^ 0. 

This is the unconstrained case discussed in Section 4.3, that is, all the com- 
mutation relations (4.16) are automatically satisfied. We can diagonahze to 
give 



W, 



(4) 



/ a; \ 



A'2 

A^j 




where 

(a;, A'2, A^) e {(1, 1, 1), (1, 1, 0), (1, 0, 0), (0, 0, 0), (1, 1, -1), (1, -1, 0)} , 

so there are six distinct cases. The exact form of the transformation is unim- 
portant, but the (1,1,0) extension can be mapped to Case 2 (the trans- 
formation is complex), (1,0,0) can be mapped to Case 3a, and (1,-1,0) 
can be mapped to Case 2. Finally the (1,1,1) extension can be mapped 
to the (1, 1, —1) case by a complex transformation. After transforming that 
(1, 1, — 1) case, we are left with 



W, 



(4) 



V 



^ 
0000 
0000 
0000 



0010 
0100 
1000 
0000 



These will be called Cases la and lb. 
Case 2 ^1 = 0, ^2 = 1- 
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The commutation relations (4.16) reduce to W^^^ = W4^^ = W4^^ = 0, and we 
have 



W, 



(4) 



( W^^^ W^^ \ 

w^^w^'' 







We can remove Wa^^ because it is a coboundary (in this case a multiple 
of ^^(3)). We can also rescale appropriately to obtain four possible exten- 
sions: VF(4) = 0, and 



(4) 



^ 1 ^ 





\ y 



1 ^ 

10 



\^ y 



1 ^ 
0-100 







Again, the form of the transformation is unimportant, but it turns out that 
the first of the above extensions can be mapped to Case 3c, and the second 
and third to Case 3b. This last transformation is complex. Thus there is only 
one independent possibility, the trivial extension Vl^(4) = 0. 

Case 3 = 1,^2 = 0. 

We can remove W/^^ using a coordinate transformation. Prom the commuta- 
tion requirement (4.16) we obtain W^'^ — W^'^ 
and 



0. We are left with Wi 



(3) 



/ W^^ Wi"^ \ 



W^(4) 
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W^^ 1^4^^ 



'4 









Using the fact that elements of the form (0, 0:25 0, 04) are an Abelian ideal of 
this bracket, we find that VF4'^'^Vl/4^^ = 0. Using an upper-triangular trans- 
formation we can also make Wi^^Wi"^ — 0. After suitable rescalings we find 
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there are five cases. One of these, 

\ 



) 

0. We are thus left with four 

\ / 1 ^ 

10 

' 10 

y \ ^ 

We will refer to these four extensions as Cases 3a-d, respectively (Case 3a is 
the trivial extension). 

Case 4 ^1 = 1, ^2 = 1- 

The elements W4^^ and W4^^ are coboundaries that can be removed by a 
coordinate transformation. Prom (4.16) we have W^^^ = W^^'^ = 0, W^^^ = 
W4^^ -. 03, so that 

^3 ^ 
^3 
^3 
^ 

For ^3 = 1 we have the Leibniz extension. The two cases will be referred to as 
Case 4a for 9^ = and 4b for 63 = 1. 

Table 1 summarizes the results. There are are total of nine independent n = 4 
extensions, four of which are trivial (Vr(4) =0). As noted in Section 4.5 only 
the Leibniz extension. Case 4b, has nonvanishing for all 1 < i < n. 

The surprising fact is that even to order four the normal forms of the extensions 
involve no free parameters: all entries in the coefficients of the bracket are 
either zero or one. There is no obvious reason this should hold true if we 
try to classify extensions of order n > 4. It would be interesting to find out, 
but the classification scheme used in this paper becomes prohibitive at such 



Wt 



(4) 



V 



1 

1 
000 
000 



may be mapped to Case 4 (below) with 63 
cases: the trivial extensions, 14^(4) = 0, and 







10 





/ 



1 



1 
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Table 1 

Enumeration of the independent extensions up to ra = 4. We have = for 

all the cases, and we have left out a row and a column of zeros at the end of each 
matrix. We have also omitted cases l-4a, for which W^(4) = 0. 



Case W, 



(2) 



Wi 



(3) 



W, 



(4) 



(0) 




^ 1 \ 

1 

1 



(0) 



1 

1 




1 

1 



V 





1 



V 



1 

10 



1 \ 

1 



1 



1 

1 
1 



high order. The problem is that some of the transformations used to relate 
extensions cannot be systematically derived and were obtained by educated 
guessing. 



5 Casimir Invariants for Extensions 



In this section we will use the bracket extensions of Section 4 to make Lie- 
Poisson brackets, following the prescription of Section 2. In Section 5.1 we 
write down the general form of the Casimir condition (the condition under 
which a functional is a Casimir invariant) for a general class of inner brackets. 
Then in Section 5.2 we see how the Casimirs separate for a direct sum of 
algebras, the case discussed in Section 4.1. Section 5.3 discusses the particular 
properties of Casimirs of solvable extensions. In Section 5.4 we give a general 
solution to the Casimir problem and introduce the concept of coextension. 
Finally, in Section 5.5 we work out the Casimir invariants for some specific 
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examples, including CRMHD and the Leibniz extension. 



5.1 Casimir Condition 



A generalized Casimir invariant (or Casimir for short) is a function C : g* 
for which 



for all F : q* 



{F,C} = 0, 
Using (2.1) and (2.4), we can write this as 



6F SC 



6C 



it 



5F_ 



Since this vanishes for all F we conclude 



5C 



t 



(5.1) 



To figure out the coadjoint bracket corresponding to (2.11), we write 
which after using the coadjoint bracket of g becomes 



SO that 



We can now write the Casimir condition (5.1) for the bracket by extension as 



6C_ 



0, u — 0, . . . ,n. 



(5.2) 



We now speciahze the bracket to the case of most interested to us, where 
the inner bracket is of canonical form (2.6). As we saw in Section 2, this is 
the bracket for 2-D fluid flows. The construction we give here has a finite- 
dimensional analogue, where one uses the Cartan-Killing form to map vectors 
to covectors, but we will not pursue this here (see Thiffeault [55]). Further, 
we assume that the form of the Casimir invariants is 



Jn 



(5.3) 



36 



and thus, since C does not contain derivatives of functional derivatives of C 
can be written as ordinary partial derivatives of C. We can then rewrite (5.2) 
as 



0, 



1/ 



0,...,n. 



(5.4) 



In the canonical case where the inner bracket is like (2.6) the C'^ are 
independent and antisymmetric in A and a. Thus a necessary and sufficient 
condition for the Casimir condition to be satisfied is 



for A, cr, = 0, n. Sometimes we shall abbreviate this as 



(5.5) 



(5.6) 



that is, any subscript on C following a comma indicates differentiation with 
respect to Equation (5.6) is trivially satisfied when C is a hnear function 
of the ^'s. That solution usually follows from special cases of more general 
solutions, and we shall only mention it in Section 5.4.2 where it is the only 
solution. 



An important result is immediate from (5.6) for a semidirect extension. When- 
ever the extension is semidirect we shall label the variables . . . , be- 
cause the subset ■ ■ ■ then forms a solvable subalgebra (see Section 4.4 
for terminology). For a semidirect extension, W^^^ is the identity matrix, and 
thus (5.6) gives 



Act ctA ) 

which is satisfied because we can interchange the order of differentiation. 
Hence, u = does not lead to any conditions on the Casimir. However, the 
variables /i, A, a still take values from to n in (5.6). 



5.2 Direct Sum 



For the direct sum we found in Section 4.1 that if we look at the 3-tensor W 

as a cube, then it "blocks out" into smaller cubes, or subblocks, along its main 
diagonal, each subblock representing a subalgebra. We denote each subblock 
of Wx^"^ by Wix'^^ 1 i = 1, . . . , r, where r is the number of subblocks. We can 
rewrite (2.1) as 
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{A,B} = Y.Ui,w,r 



i=i 

r 



5A SB 



where i labels the different subblocks and the greek indices run over the size 
of the ith subblock. Each of the subbrackets { , } . depends on different fields. 
In particular, if the functional C is a Casimir, then, for any functional F 

{F,C} = j2{F,C}, = =^ {F,C}^ = 0, t = l,...,r. 

i=l 

The solution for this is 

= C,[^i] + ■■■ + Cr[^r] , where {F , Q}. = 0, i = 1, . . . , r , 

that is, the Casimir is just the sum of the Casimir for each subbracket. Hence, 
the question of finding the Casimirs can be treated separately for each compo- 
nent of the direct sum. We thus assume wc arc working on a single degenerate 
subblock, as we did for the classification in Section 4, and henceforth we drop 
the subscript i. 

There is a complication when a single (degenerate) subblock has more that 
one simultaneous eigenvector. By this we mean k vectors u^"'\ a — 1, . . . ,k, 
such that 

Note that lower-triangular matrices always have at least the simultaneous 
eigenvector = 5^". Let r/^"^ := u'f^^'^, and consider a form C{r]^^\ . . . ,7]'^^^) 
for the Casimir. Then 



a2p k pap 



.AW ^ 



a,b=l 



Because the eigenvalue A*^*^) docs not depend on a (the block was assumed to 
have degenerate eigenvalues), the above expression is symmetric in A and a. 
Hence, the Casimir condition (5.5) is satisfied. 

The reason this is introduced here is that if a degenerate block splits into a 
direct sum, then it will have several simultaneous eigenvectors. The Casimir 
invariants C^"-'>{ri^"'^) and C'-^-' (77 '•''•') corresponding to each eigenvector, instead 
of adding as CW(r;(°)) +C^^\r]^^^), will combine into one function, C{ri^°'\ri^^^), 
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a more general functional dependence. However, these situations with more 
than one eigenvector are not limited to direct sums. For instance, they occur 
in semidirect sums. In Section 6 we will see examples of both cases. 



5.3 Local Casimirs for Solvable Extensions 



In the solvable case, when all the W^^'^'s are lower-triangular with vanish- 
ing eigenvalues, a special situation occurs. If we consider the Casimir condi- 
tion (5.4), we notice that derivatives with respect to ^" do not occur at all, 
since W^"'^ = 0. Hence the functional 

= / r (x') 5(x - x') d'x' = 
Jci 

is conserved. The variable .^"(x) is locally conserved. It cannot have any dy- 
namics associated with it. This holds true for any other simultaneous null 
eigenvectors the extension happens to have, but for the solvable case is 
always such a vector (provided the matrices have been put in lower-triangular 
form, of course). 

Hence there are at most n — 1 dynamical variables in an order n solvable 
extension. An interesting special case occurs when the only nonvanishing W(^n) 
is ior fjL — n. Then the Lie-Poisson bracket is 



{F,G} 



n-l 
/i,i/=l 



fe 



SF 



SG 



^{''(x) ' 5e(x) 



where ^"(x) is some function of our choosing. This bracket is not what we 
would normally call Lie-Poisson because ^"(x) is not dynamical. It gives equa- 
tions of motion of the form 



which can be used to model, for example, advection of scalars in a specified 
flow given by ^"(x). This bracket occurs naturally when a Lie-Poisson bracket 
is linearized [26,34]. 



5.4 Solution of the Casimir Problem 



We now proceed to find the solution to (5.4). We assume that all the W^^'K 
jji = 0, . . . , n, are in lower-triangular form, and that the matrix W^^^ is the 
identity matrix. Although this is the semidirect form of the extension, we will 
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see that we can also recover the Casimir invariants of the solvable part. We 
assume z/ > in (5.4), since z/ = does not lead to a condition on the Casimir 
(Section 5.1). Therefore Wx"'^ = 0. Thus, we separate the Casimir condition 
into a part involving indices ranging from 0, . . . , n — 1 and a part that involves 
only n. The condition 

n 

IX,(T,\=0 



becomes 



E 

A=0 



/ -I 

n—1 

E ^rc,. 

, fj,,a=0 



n—1 



e\r] + E^A'^'^c,^n[e\r 

//=0 



0, 



where we have used Wx""" = to limit the sum on /x. Separating the sum on A 
gives 



n-l 

E 

A=0 



r n-l 



n-l 



E ^a'^t, ,a[e,e] + i: wx^'c c 

, //,(T=0 /i=0 

n—1 n—1 
/i,cr=0 /i=0 



The last sum vanishes because [ ^" , ] = 0. Now we separate the condition 
into semisimple and solvable parts, 



n-l 

E 

M=i 



n-l n—1 
, \,a=0 (7=0 



n-l 



n-l 



A,o-=0 
n-l 



(T=0 

n-l 

- E On [ r , r ] + E Wn'^'c 0, [ r , r ] = 



cr=0 



a=0 



Using Wo-^'^ = ^o-'^j we can separate the conditions into a part for u = n and 
one for < 1/ < n. For u = n, the only term that survives is the last sum 

n-l 

ECoa[r,r] = o. 

(7=0 

Since the commutators are independent, we have the conditions, 

Coa = 0, (7 = 0,..., n-l. (5.7) 

and ioT < u < n, 
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n-1 

E 



n— 1 n— 1 

E ^a'^'-c [ e\ r ] - E ^-'^'c, ^„ [ c , r 

,A,cr=l cr=l 

n-1 



(T=l 



-Con[r,r] = o, 



where we have used (5.7). Using independence of the inner brackets gives 



(5.8) 
(5.9) 



for < ajXjU, 11 < n. From now on in this section repeated indices are 
summed, and all greek indices run from 1 to n — 1 unless otherwise noted. 
We have written a tilde over the 1^'s to stress the fact that the indices run 
from 1 to n — 1, so that the W represent a solvable order (n — 1) subexten- 
sion of W. This subextension does not include We have also made the 

definition 



(5.10) 



Equation (5.8) is a Casimir condition: it says that C is also a Casimir of W. We 
now proceed to solve (5.9) for the case where g is nonsingular. In Section 5.4.2 
we will solve the singular g case. We will see that in both cases (5.8) follows 
from (5.9). 



5.4-1 Nonsingular g 

The simplest case occurs when g has an inverse, which we will call g^^. 
Then (5.9) has the solution 

where 

A'^^,:=g,,W„'"'. (5.12) 

It is easily verified that Aijf^ = A^^, as required by the symmetry of the left- 
hand side of (5.11). 

In (5.11), it is clear that the nth variable is "special"; this suggests that we 
try the following form for the Casimir: 

c{e, e\ • • • , D = E ^^'^ (e", C-') MC), (5.13) 

i>0 
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where / is arbitrary and fi is the ith derivative of / with respect to its argu- 
ment. One immediate advantage of this form is that (5.8) follows from (5.9). 
Indeed, taking a derivative of (5.9) with respect to inserting (5.13), and 
equating derivatives of / leads to 

where we have used (5.7). Since the left-hand side is symmetric in A and a 
then so is the right-hand side, and (5.8) is satisfied. 

Now, inserting the form of the Casimir (5.13) into the solution (5.11), we can 
equate derivatives of / to obtain, for r, a = 1, . . . , n — 1, 



V% = Q, r,a=l,...,n-l; (5.14) 
V%^A>^,,V^:-'^ + g,„V%-'\ i>l. (5.15) 



The first condition, together with (5.7), says that V^^'' is linear in . . . 
There are no other conditions on V^^\ so we can obtain n independent solu- 
tions by choosing 

= jy = 0,...,n-l. (5.16) 

The equation for D'^^^'^ is 

"^"^ (5.17) 

Thus 'D(^)^ is a quadratic polynomial (the arbitrary linear part does not yield 
an independent Casimir, so we set it to zero). Note that D^^)^ does not depend 
on since r.a = 1, . . . ,n ~ 1. Hence, for z > 1 we can drop the term 
in (5.15). Taking derivatives of (5.15), we obtain 

-r){i)i' _ Af^i /1M2 . . . /l^f'-i) 7^(1)^^ f'f; IS"! 

^,rir2...T(i+i) - ^nra ^Mi^3 ^M(«~2)^« M(i-i)r(i+i) ■ {O.IQ) 

We know the series will terminate because the W^'^\ and hence the are 
nilpotent. The solution to (5.18) is 

pW^ ^ 1 cn^T2 . . . er(i+i) -^^ -j^gx 

(i-M)! '"i^2-T(i+i)'; , ^ ^ -l, ^--l^; 

where the constants D are defined by 

In summary, the pW's of (5.13) are given by (5.16), (5.17), and (5.19). 
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Because the left-hand side of (5.18) is symmetric in all its indices, we require 



A'' — Af^ A- 

T<7 IXA TA , 



i > 1. 



(5.21) 



This is automatically satisfied for the nonsingular g case [55] . Comparing this 
to (2.13), we see that the A's satisfy all the properties of an extension, except 
with the dual indices. Thus we call the A's the coextension of W with respect 
to g. Essentially g serves the role of a metric that allows us to raise and lower 
indices. 

For a solvable extension we simply restrict u > and the above treatment 
still holds. We conclude that the Casimirs of the solvable part of a semidirect 
extension are Casimirs of the full extension. We have also shown, for the case 
of nonsingular g, that the number of independent Casimirs is equal to the 
order of the extension. 

5.4-2 Singular g 

In general, g is singular and thus has no inverse. However, it always has a 
(symmetric and unique) pseudoinverse 51^^^ such that 



The pseudoinverse is also known as the strong generalized inverse or the 
Moore-Penrose inverse [56]. It follows from (5.22) and (5.23) that the ma- 
trix operator 



projects onto the range of g. The system (5.9) only has a solution if the 
following solvability condition is satisfied: 



that is, the right-hand side of (5.9) must live in the range of g. 

If C On 7^ 0, the quantity W^'^^C^f^n + ^'^o-C.on has rank equal to n, because 
the quantity W^"^ C_ ^„ is lower-triangular (it is a linear combination of lower- 
triangular matrices). Hence the projection operator must also have rank n. But 
then this implies that g has rank n and so is nonsingular, which contradicts 
the hypothesis of this section. Hence, C on = for the singular g case, which 
together with (5.7) means that a Casimir that depends on can only be of 
the form C = /(C'')- However, since is not an eigenvector of the W^^^^s, the 
only possibihty is C = the trivial linear case mentioned in Section 5.1. 




OtJ-cr g gru 9fJ.i^ ) 



(5.22) 
(5.23) 




(5.24) 



43 



The solvability condition (5.24) can thus be rewritten as 

(P^ W,-'^ - H^/^) C ^„ = 0. (5.25) 

An obvious choice would be to require P'^j. Wa-^^ = W^f'^^, but this is too 
strong. We will derive a weaker requirement shortly. 

By an argument similar to that of Section 5.4.1, we now assume C is of the 
form 

C{i\ C) = T.1^^He, C~') MC), (5.26) 

i>0 

where again fi is the ith derivative of / with respect to its argument. As in 
Section 5.4.1, we only need to show (5.9), and (5.8) will follow. The number 
of independent solutions of (5.9) is equal of the rank of g. The choice 

V^^> = p-^^p^ i/ = l,...,n-l, (5.27) 

provides the right number of solutions because the rank of P is equal to the 
rank of g. It also properly specializes to (5.16) when g is nonsingular, for 
then P% = 5%. 

The solvability condition (5.25) with this form for the Casimir becomes 

{P\ W/" - W^'"") = 0, i > 0. (5.28) 
For i = the condition can be shown to simplify to 

P\ W/^ = W^""^ P\, 
or to the equivalent matrix form 

PW^,)^W^,)P, (5.29) 

since P is symmetric [56]. 
Equation (5.9) becomes 

g^^V<;'>^0^ 

If (5.25) is satisfied, we know this has a solution given by 

where £ is arbitrary, and (5^^ — 9\p 9^^) projects onto the null space of g. 
The left-hand side is symmetric in A and cr, but not the right-hand side. 
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We can symmetrize the right-hand side by an appropriate choice of the null 
eigenvector, 



in which case 



where 



- 9ap Wr + 9xp W/'' - gxp Qa. g'^ W;^^ , (5.30) 

which is symmetric in A and a. Equation (5.30) also reduces to (5.12) when g 
is nonsingular, for then the null eigenvector vanishes. The full solution is thus 
given in the same manner as (5.18) by 

'^^^'^ - rr^ • • • r^'-' , ^ > o, (5.31) 

where the constants D are defined by 
and P^o) is given by (5.27). 

The A's must still satisfy the coextension condition (5.21). Unlike the nonsin- 
gular case this condition does not follow directly and is an extra requirement 
in addition to the solvabihty condition (5.28). Note that only the i — case, 
Eq. (5.29), needs to be satisfied, for then (5.28) follows. Both these condi- 
tions are coordinate-dependent, and this is a drawback. Nevertheless, we have 
found in obtaining the Casimir invariants for the low-order brackets that if 
these conditions are not satisfied, then the extension is a direct sum and the 
Casimirs can be found by the method of Section 5.2. However, this has not 
been proven rigorously. 



5. 5 Examples 



We now illustrate the methods developed for finding Casimirs with a few 

examples. First we treat our prototypical case of CRMHD, and give a phys- 
ical interpretation of invariants. Then, we derive the Casimir invariants for 
Leibniz extensions of arbitrary order. Finally, we give an example involving a 
singular g. 
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5.5.1 Compressible Reduced MHD 



The W tensors representing the bracket for CRMHD (see Section 2.2.4) were 
given in Section 2.3.2. We have n = 3, so from (5.10) we get 



[ -A| 






-A j 






In this case, the coextension is trivial: all three matrices A^^'' defined by (5.12) 
vanish. Using (5.13) and (5.16), with v — 1 and 2, the Casimirs for the solvable 
part are 

and the Casimir associated with the eigenvector is 

Since g is nonsingular we also get another Casimir from the semidirect sum 
part, 

c° = ef{e) -jeene) = a;/(^) - ^vvn^j). 

The physical interpretation of the invariant is given in Morrison [25] and 
Thiffeault and Morrison [21]. This invariant implies the preservation of con- 
tours of V', so that the value %1)q on a contour labels that contour for all times. 
This is a consequence of the lack of dissipation and the divergence-free nature 
of the velocity. Substituting C^(V') = '4>^ we also see that all the moments of 
the magnetic flux are conserved. By choosing C^{ip) = Q{ijj{'x) — cuq), a heavy- 
side function, and inserting into (5.3), it follows that the area inside of any 
ip-contouY is conserved. 

To understand the Casimirs and C^, we also let g{tlj) — — ipo) in C^. 
In this case we have 

C'^[v;tjj]= / V g{ijj) d^x = / f(x)d^a;, 

where represents the (not necessarily connected) region of Q enclosed by 
the contour — ipQ and ^^'o is its boundary. By the interpretation we gave 
of C'^, the contour moves with the fluid. So the total value of v inside of 
a ■(/'-contour is conserved by the flow. The same is true of the pressure p. (See 
Thiffeault and Morrison [21] for an interpretation of these invariants in terms 
of relabeling symmetries, and a comparison with the rigid body.) 
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The total pressure and parallel velocity inside of any ■^/'-contour arc preserved. 
To understand C^, we use the fact that lu — and integrate by parts to 
obtain 

C'iu, v,p, V'] = - (v</. • VV' + ^) fW d'x. 

The quantity in parentheses is thus invariant inside of any ^/'-contour. It can 
be shown that this is a remnant of the conservation by the full MHD model 
of the cross helicity, 

[ v-Bd^x, 

Jn 

at second order in the inverse aspect ratio, while the conservation of [v ; ip] 
is a consequence of preservation of this quantity at first order. Here B is the 
magnetic field. The quantities C^[ip] and C^\p;ip] they are, respectively, the 
first and second order remnants of the preservation of helicity, 

[ A-Bd^x, 
Jn 

where A is the magnetic vector potential. 



5.5.2 Leibniz Extension 

We first treat the nilpotent case. The Leibniz extension of Section 4.5 can be 
characterized by 

W^a^" = 5a^+% fi,u,X = l,...,n, (4.13) 

where the tensor 6 is the ordinary Kronecker delta. Upon restricting the indices 
to run from 1 to n — 1 (the tilde notation of Section 5.4), we have 

The matrix g^s nonsingular with inverse equal to itself: gi^i^ = ^/j+v The co- 
extension of W is thus 

n— 1 n— 1 

v=l i/=l 

Equation (5.20) becomes 



-^TiT2...r(j+i) ^TiT2 ^HiTs ^A'(i-2)-ri 1) -'"(j+i) 
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which, as required, is symmetric under interchage of the Tj. Using (5.13), 
(5.16), (5.17), and (5.19) we obtain the n — 1 Casimir invariants 

c'ie, • • • ' = g (TTi)! '^-+-+•••+^(^.1) • • • ^^-33) 

for z/ = l,...,n — 1. The superscript u on f indicates that the arbitrary func- 
tion is different for each Casimir, and recaU the subscript i denotes the ith 
derivative with respect to The nth invariant is simply C'^(^") = /"(C"), 
corresponding to the null eigenvector in the system. Thus there are n inde- 
pendent Casimirs, as stated in Section 5.4.1. 

For the Leibniz semidirect sum case, since g is nonsingular, there will be an 
extra Casimir given by (5.33) with u — 0, and the sums run from to n — 1. 
This is the same form as the u — 1 Casimir of the order (n -|- 1) nilpotent 
extension. 

For the ith term in (5.33), the maximal value of any tj is achived when all 
but one (say, ti) of the tj are equal to n — 1, their maximum value. In this 
case we have 

n + T2 H \- Tj+i = n + i{n - 1) = v + in, 

so that Ti — i + u. Hence, the iih term depends only on (^'^'^', . . . , ^"), and 
the uth Casimir depends on (^'', . . . , ^"). Also, 

max (ti H h n+i) = {i + l){n - 1) = v + in, 

which leads to maxi — n — u — 1. Thus the sum (5.33) terminates, as claimed 
in Section 5.4.1. We rewrite (5.33) in the more complete form 

C'iC, ■ ■ ■ , D = E Ti s'^^n r • • • e fLiin 
k=l ■ 

ior u — 0, . . . ,n. Table 2 gives the u — 1 Casimirs up to order n — 5. 

5.5.3 Singular g 
Now consider the n 



4 extension from Section 4.6.4, Case 3c. We have 



'^100^ 





9 = 



'^001^ 



1 
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Table 2 

Casimir invariants for Leibniz extensions up to order n = 5 (z^ = 1). The primes 
denote derivatives. 



n Invariant 



1 fie) 

2 ene) 



3 ef{e)+ue?f{e) 



4 ef{e)+eene)+wme) 



with W(i) = W(3) 
operator is 



0. The pseudoinverse of gi is gi ^ = g and the projection 



^10 0^ 




The solvabihty condition (5.29) is obviously satisfied. We build the coextension 
given by (5.30), which in matrix form is 



= g-^ + {W^"^ g-y - g-^ g g-\ 



to obtain 



^ 



1 
1 



^(2) ^ ^(3) ^ Q_ 



These are symmetric and obviously satisfy (5.21), so we have a good coexten- 
sion. Using (5.26), (5.27), (5.31), and (5.32) we can write, ior u — 1 and 3, 



This extension has two null eigenvectors, so from Section 5.2 we also have 
the Casimir /i(C^,C'^). The functions /, g, and h are arbitrary, and the prime 
denotes differentiation with respect to argument. 
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6 Casimir Invariants for Low-order Extensions 

Using the techniques developed so far, we now find the Casimir invariants for 
the low-order extensions classified in Section 4.6. We first find the Casimir 
invariants for the solvable extensions, since these are also invariants for the 
semidirect sum case. Then, we obtain the extra Casimir invariants for the 
semidirect case, when they exist. 

6.1 Solvable Extensions 

Now we look for the Casimirs of solvable extensions. As mentioned in Sec- 
tion 5.3, the Casimirs associated with null eigenvectors (the only kind of 
eigenvector for solvable extensions) are actually conserved locally. We shall 
still write them in the form C = fd""), where C is as in (5.3), so they have 
the correct form as invariants for the semidirect case of Section 6.2 (for which 
they are no longer locally conserved). 

6.1.1 n=l 

Since the bracket is Abelian, any function C = C(^^) is a Casimir. 

6.1.2 n^2 

For the Abelian case we have C = C{^^, ^^). The only other case is the Casimir 
of the Leibniz extension, 

c(e\a=eV(a+^(a- 

6.1.3 n=3 

As shown in Section 4.6.3, there are four cases. Case 1 is the Abelian case, 
for which any function C — C(^^, is a Casimir. Case 2 is essentially the 

solvable part of the CRMHD bracket, which we treated in Section 5.5.1. Case 3 
is a direct sum of the Leibniz extension for n = 2, which has the bracket 

[(ai,a2) , {(^1,(^2)] = (0, [ai ,/3i]), 

with the Abelian algebra [«3,/93] = 0. Hence, the Casimir invariant is the 
same as for the n = 2 Leibniz extension with the extra dependence of the 
arbitrary function (see Section 5.2). Finally, Case 4 is the Leibniz Casimir. 
These results are summarized in Table 3. 
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Table 3 

Casimir invariants for solvable extensions of order n = 3. 



Case Invariant 

2 efie)+e9ie)+h{e) 

4 efie) + u^rne) + egie) + ne) 

Cases 1 and 3 are trivial extensions, that is, the cocycle appended to the n — 2 
case vanishes. The procedure of then adding dependence to the arbitrary 
function works in general. 

6.14 n^4 

As shown in Section 4.6.4, there are nine cases to consider. We shall proceed 
out of order, to group together similar Casimir invariants. 

Cases la, 2, 3a, and 4a arc trivial extensions, and as mentioned in Section 6.1.3 
they involve only addition of dependence to their n = 3 equivalents. Case 3b 
is a direct sum of two n = 2 Leibniz extensions, so the Casimirs add. 

Case 3c is the semidirect sum of the n = 2 Leibniz extension with an Abelian 
algebra defined by [ (eta, , (/Ss, j34) ] = (0, 0), with action given by 

P(ai,a2)(/53,/34) = (0, [ «! , /^g ] ) . 

The Casimir invariants for this extension were derived in Section 5.5.3. 

Case 3d has a nonsingular g, so the techniques of Section 5.4.1 can be apphed 
directly. 

Finally, Case 4b is the n — A Leibniz extension, the Casimir invariants of which 
were derived in Section 5.5.2. The invariants are all summarized in Table 4. 

6.2 Semidirect Extensions 

Now that we have derived the Casimir invariants for solvable extensions, we 
look at extensions involving the semidirect sum of an algebra with these solv- 
able extensions. We label the new variable (the one which acts on the solvable 
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Table 4 

Casimir invariants for solvable extensions of order n = 4. 



Case 


Invariant 


la 




lb 


^ /(r) + ^ giC) + ^ HC) + 


2 




3a 




3b 




3c 




3d 




4a 




4b 









part) by In Section 5.4.1 we showed that the Casimirs of the solvable part 
were also Casimirs of the full extension. We also concluded that a necessary 
condition for obtaining a new Casimir (other than the linear case = 
from the semidirect sum was that dct W(^n) 0- We go through the solvable 
cases and determine the Casimirs associated with the semidirect extension, if 
any exist. 

6.2.1 n^l 

There is only one solvable extension, so upon appending a semidirect part we 
have 



i. 








1 \ 












\^ 


} 


V 1 


y 



Since det W(i) ^ 0, we expect another Casimir. In fact this extension is of the 
semidirect Leibniz type and has the same Casimir form as the n = 2 solvable 
Leibniz (Section 5.5.2) extension. Thus, the new Casimir is just ^°/(^^). 
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Table 5 

Casimir invariants for semidirect extensions of order n = 5. These extensions also 
possess the corresponding Casimir invariants in Table 4. 



Case Invariant 

3d em') + [ee + m') ne) + uem^) 

4b ffie) + {ee + him ne) + heierne) + hiem^') 

6.2.2 n=2 

Of the two possible extensions only the Leibniz one satisfies det W[2) 7^ 0. The 
Casimir is thus 

6.2.3 n^3 

Cases 2 and 4 have a nonsingular 1^(3) . The Casimir for Case 2 is 

Csd=eV(a+m'(a, 

and for Case 4 it is of the Leibniz form 

6.2.4 ri^4 

Cases lb, 3d, and 4b have a nonsingular W(^4y The Casimirs are shown in 
Table 5. 



7 Discussion 

Using the tools of Lie algebra cohomology, we have classified low-order exten- 
sions. We found that there were only a few normal forms for the extensions, 

and that they involved no free parameters. This is not expected to carry over to 
higher orders (n > 4). The classification includes the Leibniz extension, which 
is the maximal extension. One of the normal forms is the bracket appropriate 
to compressible reduced MHD [17,48]. 
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We then developed techniques for finding the Casimir invariants of Lie-Poisson 
brackets formed from Lie algebra extensions. We introduced the concept of 
coextension, which allows one to explicitly write down the solution of the 
Casimirs. The coextension for the Leibniz extension can be found for arbitrary 
order, so that we were able obtain the corresponding Casimirs in general. 

It would be interesting to generalize the classification scheme presented here 
to a completely general form of extension bracket [5,10]. Certainly the type of 
coordinate transformations allowed would be more limited, and perhaps one 
cannot go any futher than cohomology theory allows. 

The interpretation of the Casimir invariants can be pushed further, both in 
a mathematical and a physical sense. Mathematically, a precise geometrical 
relation between cocycles and the form of the Casimirs could be formulated. 
The cocycle and Casimirs should yield information about the holonomy of the 
system. For this one must study extensions in the framework of their principal 
bundle description [29]. Physically we would like to attach a more precise 
physical meaning to these conserved quantities. The invariants associated with 
simultaneous eigenvectors can be regarded as constraining the associated field 
variable to move with the fluid elements [25]. The field variable can also be 
interpreted as partially labeling a fluid element. Some attempt has been made 
at formulating the Casimir invariants of brackets in such a manner [6,21], and 
an interpretation of cocycles in the context of dynamical accessiblity has been 
offered [55]. 

Sufficient conditions for stability can be obtained via the energy-Casimir 
method [22,23,25], or the related technique of dynamical accessibility [26,57]. 
In both these case, we can make use of the coextension to derive the stability 
conditions for Lie-Poisson bracket extensions and a large class of Hamiltoni- 
ans [55]. 
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A Proof of = I 



Out goal is to demonstrate that through a series of lower-triangular coor- 
dinate transformations we can make l^^-*^^ (which has an n-fold degenerate 
eigenvalue equal to unity) equal to the identity matrix, while preserving the 
lower-triangular nilpotent form of 

We first show that we can always make a series of coordinate transformations 
to make W\^^ = 5x^. First note that if the coordinate transformation M is 
of the form M = I + L, where I is the identity and L is lower-triangular 
nilpotent, then W^^^ = M^^ W^^^ M still has eigenvalue 1, and for /j, > 1 
the H^(^) = M-^ W^^'^ M are still nilpotent. 

For A > 1 we have 

W^'' = Wx'' + Wy}^ L,' = Wx'' + E Wx'- L,' + Lx\ (A.l) 

where we used VFa^^ — 1- Owing to the triangular structure of the set of 
equations (A.l) we can always solve for the Lx^ to make Wx^^ vanish. This 
proves the first part. 

We now show by induction that if Wx^^ — Sx^, as proved above, then W'--^^ is 
the identity matrix. For A = 1 the result is trivial. Assume that W^^" = 6^'^, 
for n < X. Setting two of the free indices to one, Eq. (2.13) can be written 

= Wx'^ S^^ = Wx'" . 

Since W^^^ is lower-triangular the index fj, runs from 2 to A (since we are 
assuming A > 1): 

A 

Y.Wx'^'W^'^ = Wx'% 
and this can be rewritten, for cr < A, 

A-1 

Finally, we use the inductive hypothesis 

;U=2 
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which is vahd for o" < A. Hence, Wx^^ — Sx' and we have proved the result. 
(W^a'^^ must be equal to one since it lies on the diagonal and we have already 
assumed degeneracy of eigenvalues.) 
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